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Abstract. We investigate what current cosmological data tells us about the 
cosmological expansion rate in a model independent way. Specifically, we study if the 
expansion was decelerating at high redshifts and is accelerating now, without referring 
to any model for the energy content of the universe, nor to any specific theory of 
gravity. This differs from most studies of the expansion rate which, e.g., assumes some 
underlying parameterised model for the dark energy component of the universe. To 
accomplish this, we have devised a new method to probe the expansion rate without 
relying on such assumptions. 

Using only supernova data, we conclude that there is little doubt that the universe 
has been accelerating at late times. However, contrary to some previous claims, we can 
not determine if the universe was previously decelerating. For a variety of methods used 
for constraining the expansion history of the universe, acceleration is detected from 
supernovae alone at > 5a, regardless of the curvature of the universe. Specifically, 
using a Taylor expansion of the scale factor, acceleration today is detected at > 12a. 
If we also include the ratio of the scale of the baryon acoustic oscillations as imprinted 
in the cosmic microwave background and in the large scale distribution of galaxies, 
it is evident from the data that the expansion decelerated at high redshifts, but only 
with the assumption of a flat or negatively curved universe. 
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1. Introduction 

It is becoming generahy acknowledged that the observed redshift-distance relation from 
Type la supernovae (SNe la), together with the scales of baryon acoustic oscillations as 
observed in the distribution of galaxies on large scales (BAO) and the temperature 
anisotropics in the cosmic microwave background (CMB), implies that the current 
energy density in the universe is dominated by dark energy, here defined as a component 
with an equation of state, w = p/ p < —1/3. Current cosmological data are consistent 
with the standard - or concordance - cosmological model, where Qm = 0.3 and Qa = 0.7, 
i.e., the dark energy is explained in terms of a cosmological constant or vacuum energy 
with w = —1 [e.g., m HU 1151 HSl [211 ES]- For the concordance model, the universal 
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expansion is accelerating at redshifts lower than 2; ~ 0.7 and decelerating at higher 
redshifts. 

Most analyses which attempt to infer something about the expansion history of the 
universe rely on a specific model, e.g., a dark energy 'fluid' or a modification of general 
relativity, which has one or two parameters of interest [13]. Consequently, much of our 
knowledge of the expansion history of the universe has these parameterisations hard- 
wired into our conclusions, which may leave a large space of possible expansion histories 
unexplored. Alongside these analyses, it is therefore constructive to try to assert model 
independent statements where we can. The question we address in this paper is to what 
extent we can infer changes in the expansion rate without referring to any theory of 
gravity or model for the energy content of the universe. That is; what is the history of 
the universal expansion? 

Unfortunately, it is difficult to give a definite answer to this question since two 
important assumptions are implicit in all discussion of this kind. It is known that 
spherically symmetric void - or Hubble bubble - models may explain the anomalous 
Hubble diagram whilst always maintaining a decelerating expansion rate, at the price 
of violating the Copernican principle [SlIllElIinillSlIISlEniESlllSlllH]. Ideally 
we should test the Copernican principle in a model independent way [H HS], and 
so rule these models out. A further assumption in the standard model - insofar as 
determining the expansion dynamics is concerned - is that that the universe is smooth 
enough at small distances to be described by a perfectly homogeneous and isotropic 
model [ai5l[T2l|26l|271|28l[32l[33l[37lllT]. At best this gives a small error to all our 
considerations; at worst, many of our conclusions might be wrong. 

Since it was traditionally thought that the expansion rate would be decelerating, 
we measure acceleration with the deceleration parameter, g, defined by 
aa a 
^"-^ = -^' 

where a is the scale factor, dots denote derivatives with respect to time and the Hubble 
parameter, H, is defined as H = a/ a. Because of the sign convention, negative values for 
q correspond to acceleration. Traditionally, the goal when observing the expansion of the 
universe was to constrain two parameters, the current values of the Hubble parameter 
- or the Hubble constant - Hq, and the deceleration parameter, go- The first signs 
of an accelerated expansion came 10 years ago with the observations that distant SN 
la appear dimmer than expected in a universe with constant or decelerated expansion 
velocity [301 IMl El] . However, the acceleration was only evaluated in terms of a model 
with rim and in which go = — f^A, see Eq. Q. 

In 2002, Turner and Riess [33] studied the change in the expansion velocity - 
without referring to the energy content or theory of gravity - by using a step model 
for the deceleration parameter where g had one constant value at low and intermediate 
redshifts and another constant value at high redshifts. They demonstrated, that the SN 
la data at the time showed a strong preference for acceleration today and deceleration 
in the past, if the transition redshift was set (by hand) to z = 0.4 — 0.6. However, this 
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conclusion relied heavily on the observed magnitude of a single supernova, SN1997ff, 
the most distant SN la observed at z = 1.755. Being very distant, it is especially 
susceptible to systematic effects that may reduce its cosmological utility. One such 
effect is gravitational lensing that has been shown to brighten SN1997ff by ~ 0.15 
magnitudes [23] . 

Also, in Shapiro and Turner |13], it was shown that marginalising over the transition 
redshift, considerably relaxed the constraints on the expansion history. Specifically, 
using the so called "gold" dataset consisting of 157 SNe la [36], the authors only found 
strong evidence for acceleration at some epoch, not necessarily at z < 0.1, and that q 
was higher in the past. 

In 2004, Riess et al. [36] , the gold dataset was used to constrain a Taylor expansion 
of g(z), 

q{z) = qo + zq'{z = 0) , (2) 

where the prime denotes a derivative with respect to redshift. The result found was 
that go ^ —0.3 and g' > at 95 % confidence level (CL). The evidence that q' > 0, 
was interpreted as evidence for deceleration at higher redshifts. However, the Taylor 
expansion should only be meaningful for z < 1, up to which the parameter space [go, q'] 
still allows for acceleration at 95% CL. Also, the fit is done using data at z > 1, for 
which higher order terms in the Taylor expansion should be important. 

In, Elgar0y and Multamaki [T7], the same data was found to be consistent with 
a constant negative deceleration parameter, and in Rapetti et al. [3_lj, the SN data 
was combined with X-ray cluster gas mass fraction measurements to constrain the 
deceleration parameter, g, as well as the next order derivative of the scale factor as 
decoded in the jerk parameter, j. See also Daly et al. [H] for an alternative approach 
to constraining the acceleration history of the universe. 

In this paper, we follow the spirit of previous work, in that we strive to make 
as few assumptions as possible regarding the theory of gravity or the energy content 
of the universe when inferring the state of the expansion velocity of the universe. In 
fact, the only assumptions used in this paper, other than those mentioned, is that SNe 
la are standardisable candles and that the observed inhomogeneities in the large scale 
distribution of galaxies and the anisotropics in the temperature of the CMB reflects 
the same physical scale. One feature we discuss in particular is the role of curvature 
in determining the acceleration, as there are significant degeneracies between curvature 
and acceleration, in a similar vein which exists between curvature and the dark energy 
equation of state w [HI [21] . 

In Sec. [2| we discuss acceleration and deceleration within the standard model, as 
well as possible observational measures of the deceleration parameter. In Sec. [3| we 
present the two sources of data used in this paper, and in Sec. |4] we present a new 
method for inferring the expansion history of the universe, together with our results for 
q{z). Our results are summarised in Sec. |5} In short, we conclude that the evidence for 
late time acceleration of the universal expansion is very strong, regardless of the method 
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Figure 1. The range of behaviour available to q{z) for the ACDM model. 



used to measure q{z), whereas deceleration at high redshifts can still be avoided given 
current data. 

2. Acceleration and deceleration 

First, let us calculate what to expect for q{z) in the standard ACDM model, with matter 
density, fi^, and a cosmological constant, f^A = 1 — Vt^ ~ ^k- The acceleration equation 
is given by 

d AnG , „ , 47rG , ^ . 

- = —{e + 3p) = — (p^-2pa). (3) 



This gives 



1 nmii + zf -2{i-nm-n 



. (4) 

In the infinite future {z = —1), we have q = —1, and as z — > ocjH q = 0.5. The transition 
redshift between acceleration and deceleration is given by zt = [2(1 — rim) /^mY^^ — 1 
when = 0. For the flat concordance cosmology with Q^n = 0.3, we have Zt ~ 0.7. One 
of the main problems we have in determining acceleration as a function of redshift, is the 
wide range of behaviours q{z) can exhibit, especially if curvature is present. In Fig[l| 
we show the range of values q{z) can take on within the ACDM paradigm by varying 
the constants Qm and Qk- With more exotic dark energy models the variation can be 
much more elaborate and q{z) can be regarded as a free function to be constrained by 
observations. Then, once q{z) is determined, the dark energy equation of state, w{z), 
can in principle be determined from a first-order differential equation. However, since 
the integration constant (i.e., Hq), as well as Qk and Qm sue arbitrary, there is a three 
parameter family of dark energy models which can give rise to the same q{z). 



I In fact, this approximation is only valid after the universe became matter dominated at z ^ 3000. 
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we may derive 



+ - 1 = (1 + z) 



In 



H 



(5) 



Normalising the scale factor to be unity today, a = {1 + z) 
z = —(1 + z)H, and 

In principle, Eq. ^ can be used to obtain the value of the deceleration parameter as a 
function of redshift. Unfortunately, it is very difficult to measure the Hubble parameter, 
H{z), let alone its derivative. As an example, from Fig. [3| it can be understood why it 
is difficult to differentiate noisy SN la data in order to obtain H{z). Note however that 
with future BAO data, it may be possible to measure H{z) directly ^2\ . 
Rearranging Eq. ([s]), we can write 



H' 



l + z 



d{lnH) 

dz 



dz 



1 + q{z) 

l + z 



dz 



If the universe is accelerating between redshifts zi and Z2, we have 



In 



El 
Hi 



< In 



Z2 



l + Zi 



{I+Z2) ^ (l + ^l) 



> 



(6) 



(7) 



H2 Hi 

When we have acceleration, H{z) thus grows with redshift faster than (l + z). This is 
easy to understand since if d = H/{1 + z) is increasing with redshift, a is decreasing 
with time and d < 0, corresponding to acceleration. For deceleration, H{z) grows slower 
than (l + z). For a matter dominated universe, H = Ho^nmil + zy oc (1 + i.e., 
the expansion is decelerating. For an empty universe, H = Ho^n^il + zy oc (1 + z), 
i.e., the expansion velocity is constant. Note however that an empty universe is not the 
only possible solution for a constant expansion velocity - any universe where the total 
energy density scales as (1 + z)^ also gives a constant expansion. 
The comoving coordinate distance is given by 



dr. 



dz 
Wz] 



1 

Ho 



exp 



^ [1 + q{u)]du 



'1 



u 



dv . 



(8) 



and the luminosity distance, which is the relevant quantity for SN la observations, is 
given by 

l + z 



diiz) = ^^ sin ^/-^}kHodc{z) . (9) 

The angular diameter distance, relevant for the scale of BAO and CMB, is given by 
dA = di/ {l+zY- Although having the same constant expansion rate, we would therefore 
measure different luminosity and angular diameter distances in the flat and open non- 
accelerating cases. Determining acceleration at a given redshift is then the same as 
determining if the function [{1 + z)D' / HQ^Jl — \Vlk\D'^]^ where D{z) = dL{z)HQ/{l + z), 
is an increasing function of redshift. 



3. Data 



In the last decade, there has been a formidable progress in using cosmological data to 
constrain the expansion history and the energy content of the universe. Observations 
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include, but are not restricted to, SNe la, BAO, CMB, weak gravitational lensing and 
galaxy cluster number counts. In this paper, we make use of SN la data together with a 
combination of CMB and BAO observations that only depends on the expansion history 
of the universe, not the energy content. 



3.1. Type la supernova data 

The UnionOS data set |25] is a compilation of SNe la from, e.g., the Supernova Legacy 
Survey, ESSENCE survey and HST. After selection cuts, the data set amounts to 307 
SNe la, spanning a redshift range of < z < 1.55, analysed in a homogenous fashion 
using the spectral-template-based fit method SALT. 



3.2. Baryon acoustic oscillations and the cosmic microvawe background 

The distances measured to the CMB decoupling epoch at ~ 1090 and to the BAO 
aX z = [0.2, 0.35] depend on the physical scale of the acoustic oscillations at decoupling, 
and thus on the matter and baryon density. Specifically, the position of the first peak 
in the CMB power spectrum, which represents the angular scale of the sound horizon 
at decoupling, is given by, 

- ,Mf4i±^ (10) 

where the comoving sound horizon at recombination, 

depends on the speed of sound, c^, in the early universe. The observed scale of the BAO 
is given by rs{z^) / Dy, where the so called dilation scale, Dy, is combined from angular 
diameter and radial distances according to 

1/3 

(12) 



Dy{z) 



Since we want to infer the expansion history with minimal assumptions regarding 
the energy density, we take a conservative approach and use the ratio of the observed 
scales in the CMB and the BAO. This ratio does not depend on the physical size of 
the sound horizon at decoupling but only on the assumption that the BAO and CMB 
reflects the same physical size. Percival et al. p9] derives 
dA{z^){l + z^) 



Dv{z = 0.2) 
dA{z^){l + z^) 



19.04 ±0.58 

10.52 ±0.32. (13) 
Dv{z = 0.35) ^ ^ 

The measurements (from the 2dFGRS and SDSS, respectively, combined with 3 year 
WMAP data) are correlated with correlation coefficient p = 0.39. Using WMAP 5 year 
data instead of 3 year data gives close to identical results, when combined with the BAO 
data. 
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Figure 2. Best fit values of dA{z*)/DY{z) at z = 0.2 (red) and z ~ 0.35 (blue), 
together with their la errors, as measured from the scale imprinted in the CMB and 
BAO, shown in the parameter space for ACDM. 



The constraints given by Eq. (13), individually place very weak constraints on 
acceleration. Also, these constraints are highly degenerate with the curvature. In 
Fig. |2| we show the CMB/BAO constraints in the ilk — Qo parameter space for the 
ACDM model, with their la errors. It is clear that they require a large, negative, 
deceleration parameter to be consistent with each other at Icr CL, irrespective of the 
curvature. Note also that a flat model is ruled out at this level. At 2a CL however, 
individual constraints on go and Qk becomes very weak, and the combined constraints 
are limited by the strong degeneracy exhibited between and go- Nevertheless, the 
data unambiguosly shows that Qk > —0.1, i.e., the universe is not strongly overdosed. 



4. Methods and results 



In this section, we discuss how to determine the change of the expansion rate, or 
q{z), without assuming a model for the energy content of the universe or a specific 
gravitational theory. We are thus limited to using geometrical data as opposed to 
methods that are sensitive to the growth of structure in the universe. 

First, we discuss some difficulties that occur when trying to find q{z) by comparing 
observed magnitudes of SNe la directly. This method has the advantage that it 
can, in principle, detect a change in the expansion rate without somewhat ad-hoc 
parameterisations of q{z). Some of these difficulties have been realised by Seikel and 
Schwarz |?T] - in particular the reliance on low redshift SNe la. 

Alternatively we may parametrise the expansion with piecewise, constant 
accelerations, i.e., we assume that q{z) can vary between redshift bins, but is constant 
within the bins. This is the approach introduced by Turner and Riess [H]. We 
investigate this further below, and find that it suffers from strong degeneracies with 
curvature when including CMB data. We also investigate alternative parameterisations 
of q{z), displaying the same degeneracies. However, we show that SN la data alone. 
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show that q{z = 0) < 0, irrespective of curvature. 

We then present a method for determining acceleration, first employed in Riess 
et al. [35], that relies on the fact that = {l + z)/H{z) is an increasing function when 
the expansion is accelerating which seems to the case at low redshifts. However, since 
the method relies on differentiating noisy SN la data, results at high redshifts have very 
large uncertainties. Finally we present a new method which relies on a Taylor expansion 
around an arbitrary 'sliding' redshift. 



4.I. Inferring acceleration from m{z) 

In SN la cosmology, the validity of a given cosmological model is tested by comparing 
observed peak SN la magnitudes with the theoretical magnitudes for the given model. 



m{z) = M + 5 log 



10 



IMpc 



+ 25, (14) 



where M is the absolute SN la magnitude. The normalisation of the magnitude 
(containing, e.g., M and Hq) is usually marginalised over, and constraints are derived 
by examining the redshift evolution of the magnitudes. 

Often, one presents SN la data in the form of the difference between the observed 
peak magnitudes, and the theoretical magnitudes in an empty universe, me, where 
the difference is normalised to be zero at low redshifts. Since an empty universe is 
neither accelerating nor decelerating, it is sometimes claimed that one can infer the 
state of the universal expansion from a quick visual inspection of this difference. One 
common claim (at talks and discussions, if not in papers) is that negative values of 
this difference, m — me < 0, shows that the universe is decelerating and positive values 
that it is accelerating. Another inconsistent, but nevertheless common claim, is that 
if the difference increases with redshift, the universal expansion is accelerating. If the 
difference is decreasing, the expansion is decelerating. That is, the claim is that we can 
infer the state of the expansion velocity by studying the sign of the derivative of the 
difference with respect to redshift. Let us investigate these claims. 

For the non-accelerating case, if oc [1 + z), and 

4,n = ^ln(l + 2), (15) 
-no 

and 

aL,n = sm 

Hoy/— ilk 

where subscript n refers to non-accelerating expansion. For a flat, non-accelerating, 
universe 

d,j=^-^\n{l + z). (17) 
For an empty, non-accelerating, universe, fi^ = 1, and 



^^ln(l + z) , (16) 
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For small curvature {VLk ^ 1), defining x = y/—VLkdc, we can use a Taylor expansion in 
X and write 

Hodi 1 

sm (x) (19) 



l + Z 



k 



^ ^a; - ^ ) + Oix") (20) 



~ ffo (4 + y rfi?) . (21) 

The difference between the observed SN la magnitudes and those expected in a non- 
accelerating universe is given by 

Am = m — rrir, = In ( — ) . (22) 



In 10 \dL,n 

In Seikel and Schwarz |10], the inequality (valid in a flat universe 



" [l + q{u)]du 



dv (23) 



Ho Jq [l + z) Ho 

was considered as evidence (~ 5a CL) for some period of acceleration up to redshift 
z, i.e., g < at some redshift. Allowing for spatial curvature, the evidence becomes 
signiflcantly weaker, or 1.8cr. We note [Eq. (18)], that this is equivalent to studying 
the sign of Am, if the universe is flat. Assuming we can get rid of the dependence on 
Hq and M by normalising the difference to be zero at low redshifts, we can study the 
dependence of the curvature term by looking at the difference when subtracting empty 
and flat non-accelerating cosmologies (Fig. |3]). It is obvious that Am in fact is positive 
for z ~ 0.5, regardless of the curvature of the universe. A similar comparsion between 
SNe la at low and mid/high redshifts was recently used in Seikel and Schwarz [41] to 
provide a calibration-independent test of the accelerated expansion of the universe, the 
conclusion being that the universe has accelerated at some epoch at ~ 4cr CL. From 
Eq. ( 23 ) , it is evident that the sign of Am only tells us whether the integrated expansion 
up to z is accelerating or decelerating on average, not the state of acceleration at a given 
redshift. 

We therefore turn to the second common claim, namely that if the difference Am 
increases with redshift, the universe is accelerating at that very redshift. The derivative 
of Am with respect to redshift is given by 

(A™)' = tAt (? - (25) 



In 10 \dL dL,n, 

We note that (Am)' is independent of Ho and M. For (Am)' > 0, we have 

^>^. (26) 
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Figure 3. The differences between observed SN fa magnitudes and the expected 
magnitudes for empty and flat non-accelerating universes. The differences are 
normalised to be zero at low redshifts. Error bars correspond to 68.3% CLs. 



Assuming a flat universe, this amounts to 
r dz 

H{z) / ]^<(1 + ^)M1 + ^)- (27) 



From Eq. (25), it is evident that (Am)' does not depend only on the state of acceleration 
at the given redshift, but also on the integral of the expansion. The sign of (Am)' is 
thus not directly related to the state of acceleration. As a specific example, consider the 
concordance model with Qm = 0.3 and Qa = 0.7 in Fig. |4} The lines corresponding to 
(Am)' and the deceleration parameter does not cross the zero line at the same redshift. 
We conclude that visually judging the state of the universal acceleration from a SN la 
Hubble diagram is non-trivial. 



4-2. Piecewise constant q{z) 

We divide our redshift range in [zq = 0, zi, Z2, z^, . . .], where between redshift, Zi-i 
and Zi, we have a constant deceleration parameter gj. In that bin, we then have 
H{z) oc (1 + z)'^^^^. Since the normalisation of the Hubble parameter is marginalised 
over, we can put H{zq) = Hq = 1. For < z < zi, we then have H{z) = (1 + zy'^'^^, for 
zi< z < Z2, we have H{z) = (1 + Zi^^^'^^il + z)«2+i^ g^c. 

The simplest case of a piecewise constant q{z) would be to have a single constant 
value of the deceleration parameter, i.e., q{z) = go and H{z) = Hq{1 + z)'^°^^, at all 
times. Actually, such a simple model gives a reasonable fit to SN la data, especially 
when allowing for curvature, see left panel of Fig. |5} To fit the CMB and BAO data 
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Figure 4. (Am)' — d{Am)/dz compared to q{z) for the concordance model. The 
sohd and dotted hnes shows the derivative of the difference between the concordance 
and the flat and empty non-accelerating models, respectively. There is not a one-to-one 
correspondence between the sign of (Am)' and the state of acceleration. 

using this simple model however, requires a large amount of fine-tuning since the quality 
of the fit depends very sensitively on the exact values of Vtk and go- Also, since CMB and 
BAO prefers a close to fiat solution, the model provides a poor fit to all data combined. 

We next turn to an analysis along the lines of Turner and Riess [44j where we 
assume that the deceleration parameter q{z) has one constant value at z < Zt and one 
constant value at z > Zt- We first choose Zt = 0.7 in order to maximise our chances 
to detect a difference in the change of the expansion rate, since this is the transition 
redshift for the concordance cosmology that we know provides a good fit to the data. 
Results using SN la data are shown in the right panel of Fig. [5j Allowing for curvature 
increases the size of the contours by a factor of ~ 2. Incorporating also CMB and BAO 
data has the effect of constraining the curvature down to the level where the combined 
results are very similar to SN la results when assuming zero curvature, except for a 
tail of very low values for q{z > 0.7) if the universe has a small positive curvature 
{ilk < 0, see Fig. [o]). Note that the data does not show any evidence for deceleration 
at z > 0.7. However, the evidence for acceleration at z < 0.7 is very strong; ~ 5a for 
SN la data only and ~ 7a including CMB and BAO data or assuming a flat universe. 
However, as shown in Shapiro and Turner [13], marginalising over the transition redshift, 
considerably relaxes the constraints on the expansion history. In the left panel of Fig. [6| 
we show how the confldence contours change with the transition redshift, Zt, assuming 
a flat universe. For high Zt, results are mostly sensitive to q{z < Zt), and vice versa. 
In the right panel of Fig. [6} results when marginalising over the transition redshift 
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Figure 5. Left panel: SN la constraints on a model where the deceleration parameter 
has a constant value q{z) ~ qq. Right panel: SN la constraints on a model where the 
deceleration parameter q has one constant value at z < 0.7 and one constant value at 
z > 0.7. Allowing for curvature increases the size of the contours by a factor of ^ 2. 
The evidence for acceleration at z < 0.7 is ~ Ta and ~ 5a, respectively. Contours 
correspond to 68.3%, 95.4% and 99.7% CLs. 



in the interval 0.1 < Zt < 1.0 are shown, with and without the assumption of a flat 
universe. The upper, right quadrant, corresponding to expansion histories without any 
acceleration, can be ruled with high confidence, independent of curvature. However, all 
other possibilities are still viable, including deceleration at low redshifts and acceleration 
at higher redshifts, if the transition redshift is low enough [zt ~ 0.1). 

Next, we allow for a finer redshift resolution in q{z). Note that the qi are correlated, 
but we can decorrelate the qi estimates by (following Ref. [22]) changing the basis 
through an orthogonal matrix rotation that diagonalises the covariance matrix. This 
corresponds to applying a weight function to the qi to obtain decorrelated Qi. These are 
linear combinations of qi where the weight function quantifies the redshift dependence 
of Qi. First, we use SN la data only to constrain q{z) in bins z = [0., 0.3, 0.6, 1., 1.8], 
assuming a fiat universe. Results for q^ and the decorrelated Qi together with the 
corresponding weights are presented in Fig. [7j It is interesting to note that we again 
detect acceleration at low redshift with high confidence, but that the data do not require 
deceleration at higher redshifts. Allowing for curvature increases the error bars on q{z) 
somewhat, but does not change our qualitative results. 

Including CMB and BAO data, we can extend our analysis all the way out to 
z ~ 1100 in redshift bins z = [0., 0.5, 1., 1.8, 1100.]. We can see (Fig. g) that in the 
redshift interval 1.8 < z < 1100, we in fact do see that the expansion was decelerating, 
assuming a fiat universe. From the weight function, we can also see that the deceleration 
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Figure 6. SN la constraints the deceleration parameter q when varying the transition 
redshift z^. In the left panel, confidence contours for three discrete values of Zf in a flat 
universe is shown. In the right panel, results when marginalising over the transition 
redshift in the interval 0.1 < Zt < 1.0 are shown, with and without the assumption of 
a flat universe. Contours correspond to 68.3%, 95.4% and 99.7% CLs. 
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Figure 7. SN la constraints on q{z) in bins z — [0., 0.3, 0.6, 1., 1.8], assuming a 
flat universe. The left panel shows results for qi, and the right panel decorrelated Qi 
together with the corresponding weights. Error bars represent 95.4% CLs. 
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Figure 8. Constraints on q{z) in redshift bins z — [0., 0.5, 1., 1.8, 1100.] combining 
SN la, CMB and BAO data, assuming a flat universe. The left panel shows results 
for qi, the right panel decorrelated Qi and the corresonding weight functions. In the 
redshift interval 1.8 < z < 1100, the expansion is decelerating at high CL. Error bars 
represent 95.4% CLs. 



parameter at different redshifts are very much correlated. If we relax the assumption 
of flatness, our Monte Carlo Markov Chains fail to converge, the reason being that 
the curvature and the deceleration parameter at high redshifts exhibit very interesting 
degeneracies, that also allow for solutions with high redshift acceleration and positive 
curvature. This degeneracy is easy to understand when considering the two component 
model where the deceleration parameter q{z) has one constant value aX z < zt and one 
constant value at 2; > Zt- For Zt > 0.35, and a fixed q{z < Zt), BAO constraints on 
Dy only has a very weak (sub-percent) dependence on curvature, and the ratio between 
the scale of the BAO and CMB solely depends on the angular diameter distance to 
~ 1090, 

dAiz^){l + z*) = —^===sm ^/ -ilkHodciz^) . (28) 

In the left panel of Fig. [9} we have plotted the difference between the measured and 
theoretical values of dA^z^) for a fixed value of Dv{z = 0.2), normalised with the error on 
dA^z^:). The characteristic degeneracy structure arises because of the sinusoidal form of 
the angular distance in Eq. (28 ) and allows for models with high redshift acceleration and 
positive curvature to fit the data. In the right panel, the full confidence contours for SN 
la, CMB and BAO data are shown. It is evident that negative values for q{z > 0.7), i.e., 
corresponding to acceleration, are allowed if the universe has a small positive curvature. 
We note however, that such solutions require a considerable amount of fine tuning of 
the deceleration and curvature parameters. 



Model independent constraints on the cosmological expansion rate 15 



AdA(z.)/CT[dA(z.)] 




-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 -1-5 -1-0 -0.5 0.0 0.5 1.0 

q(z>0.7) q(z>0.7) 



Figure 9. An example of the strong degeneracy between the curvature and the 
deceleration parameter at high redshifts in the simple model where q{z) has one 
constant value at z < 0.7 and one constant value at z > 0.7. In the left panel, 
the difference between the measured value and theoretical values of c?^(z*) for a fixed 
value of Dy(z = 0.2), normalised with the measured error on is shown. In the 

right panel, full confidence contours employing SN la, CMB and BAO data are shown. 
The degeneracy is not sensitive to the exact redshift binning of q{z). The shaded areas 
represent 68.3%, 95.4% and 99.7% CLs. 



4-3. Parameterising q{z) 

There are many ways to parameterise q{z). Clearly, a piecewise constant function as 
used above only gives a limited amount of information about q{z); its average value over 
a certain redshift range. Ideally, we would prefer something which can closely mimic 
what the 'true' q{z) might be up to. As an example, let us consider 

q{a) = qo + qa{l - a) = qo + qa—— , (29) 

i. -\- z 

adopted from one of the most common parameterisation of the dark energy equation of 
state, w{z) [7]. At zero redshift, q{a = 1) = go and in the infinite past, q{a = 0) = qo+qa- 
This parameterisation appears to be reasonably flexible in the sense that performing a 
least squares fit to many random dark energy models, always provides an acceptable fit. 



In Fig. [T0| we show results for this parameterisation, both with and without curvature. 
Once again curvature degrades all constraints, though not as significantly as for other 
models for the SN la constraints on q{z). Therefore, fitting SN la data alone using 
this parameterisation, we are able to conclude that the universe is accelerating today at 
> 5cr, irrespective of curvature. 

For qa = 0, the model has q{z) = qo and is equivalent to the model presented in 
Fig. |5j Studying the contours at = in Fig. 10 , it is again obvious that this simple 
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Figure 10. Constraints on la for the parameterisation given by Eq. (29 1, assuming 



a flat universe (left) and allowing for curvature (right). The dash-dotted line is given 
by (7o + 'Za = 0, and separates regions with eternal acceleration (lower left) from regions 
with past deceleration (upper right). Using only SN la data, < Q a,i > ba CL, 
irrespective of the spatial curvature. The degenerate structures in the right panel are 



discussed in Sec. 4.2 Regions represent 68.3%, 95.4% and 99.7% CLs. 



model does not provide a good fit to the data, unless we allow for curvature and an 
extreme fine-tuning of the parameters of the model. 

4.4- The Hubble parameter 

Since it is possible, although difficult (see Fig. |3]), to obtain the Hubble parameter 
by differentiating distances, which would give us q{z) directly, we can try to compare 
= {l + z)/H{z) at different redshifts |35]. We follow the technique proposed in Wang 
and Tegmark |17] to extract the expansion history in uncorrelated redshift bins from 



SN la data. The results are shown in Fig. 11, where increasing values of (1 + z)/H{z 



corresponds to acceleration and vice versa. Unfortunately, because of the differentiation 
of sparse and noisy data, results are quite sensitive to the employed binning of the data. 



especially at high redshift. In the left panel of Fig 11 , redshifts bins of size Az = 0.3 has 



been used whereas in the right panel, Az = 0.25. It is therefore difficult to assess the 



exact CL for the trend seen in Fig. 11 of increasing (1 + z)/H{z) - corresponding 
to acceleration - at low redshifts and decreasing (1 + z)/H{z) - corresponding to 
deceleration - at high redshifts. However, independent of the bin size, acceleration 
at low redshifts seems to be inevitable whereas deceleration at higher redshifts is not. 
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Figure 11. SN la constraints on d^^ = (1 + z)/H{z). Values increasing with 
redshift correspond to acceleration and vice versa. The dotted lines correspond to a 
flat universe with f2m = 0.2 and = 0.4 from top to bottom. The normalisation 
is arbitrary and is set to agree at low redshifts. In the left panel of Fig 11 redshifts 
bins of size Az = 0.3 has been used whereas in the right panel, Az = 0.25. Vertical 
error bars correspond to 68.3% CLs. Acceleration at low redshifts is clearly detected 
whereas deceleration at higher redshifts is not. 



4-. 5. Sliding Taylor expansion 



We can Taylor expand the scale factor around the current value according to 
a{t) = ao + ao(t - to) + ^ao(t - ^o)^ + 0{t - t^f 



1 + HQ{t - t 



oj 



qoHl{t-t,f + 0{t-tof. 



Since the comoving coordinate distance is given by 

"^o dt 



we can write 



4 (to 



driz) 



a{t) 



1 



1 + go 



(30) 
(31) 

(32) 

(33) 



To first order, the distance is given by the expansion rate of the universe today, or 
Hq, and to second order by the change in the expansion rate today, or go- We now 
generalise Eq. (33) to allow for a Taylor expansion around any time, or equivalently, 
redshift according to 



a(t) ^ a(ti) + di(t - ti) + -di(t - tif 



ai 



l + Hi{t-ti 



:qiHl{t-t,] 



The comoving coordinate distance is given by 

1 



4 (to 



(to - ti) 



^(to-tl)'-(te-ti) + ^(te 



(34) 
(35) 

(36) 
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Figure 12. Fitting the deceleration parameter, q{zi) — qi, with a shding expansion 
redshift for data simulated within a flat, cosmological constant universe with ri,„ — 0.3. 
The bias between the fitted qi and the true q{z) (dotted line) caused by higher order 
terms in the Taylor expansion lies within the 68.3% CL for < 0.6. 



or, in terms of observables, 



dr. 



z 

Hi 



1 + 



zi{l + qi] 



l (l + gi) 
2(l + zi)' 



We can now fit the parameter gi, using SN la data for a sliding expansion redshift zi, 
while marginalising over Hi. In Fig. 12, we show results obtained using simulated data, 
corresponding to the UnionOS data set, within the concordance cosmology. We only 
include SNe for which |2; — 2i| < 0.75 to guarantee convergence, which is why this and 
subsequent plots only extend to z = 0.75. As can be seen in Fig. 12, the bias from 
higher order terms is within the 68.3 % CL for z < 0.6 and within 95.4 % CL all the way 
up to 2; = 0.75. 

In Fig. 13, results for the real UnionOS data set, with and without curvature 
(— 1 < < 1), is shown. Positive VLk pushes q{z) towards more positive values and 



^1 



< 0.75 to reduce bias from 



vice cersa. Again, we only include data for which | 
higher order terms in the Taylor expansion. 

Regardless of the expansion redshift, the sliding Taylor expansion is able to provide 
a very good fit to the SN la data. Note however that the results for qi at different 
redshifts are not independent, and it is therefore not possible to interpret the results as 
giving a full functional form of deceleration parameter q{z). Assuming zero curvature, 
the universe can be shown to be accelerating at 2 < 0.5 and there is weak (95.4%) 
evidence for decelaration at high redshifts. Including curvature severly degrades these 
constraints, especially the behaviour of q{z) at high redshifts. Nevertheless, using the 
sliding Taylor expansion, the evidence for acceleration at z = is > 12cr, even when 
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Figure 13. Fitting the deceleration parameter, q{zi) — qi, with a shding expansion 
redshift. Note that constraints on qi are significantly looser when allowing for curvature 
(in this case — 1 < fi^ < 1). However, the evidence for acceleration at z = is still 
> 12(7, including curvature. For a flat universe, qi{zi — 0) < at ~ 15a CL. The 
lines correspond to the flat, cosmological constant universes with fi^ = 0.2 and 0.4, 
respectively. The shaded areas correspond to 68.3%, 95.4% and 99.7% CLs. 

allowing for curvature. 

The reason for the closeness of the 2 and 3a contours in the right panel, is that 
at each redshift, the confidence contours in the [Qk, (jij-plane has a banana like shape. 
This has the effect that the resulting confidence levels for qi, which are obtained by 
projecting the banana shaped contours on the gi-axis, has a non-gaussian shape. 

5. Summary 

In this paper, we have investigated to what extent we can measure the change in the 
universal expansion rate, without making any assumptions about the energy content of 
the universe. Consequently we are limited to geometrical data as opposed to methods 
that are sensitive to the growth of structure in the universe. The data employed in this 
paper includes the redshift-distance relation of Type la SNe, as well as the the ratio 
of the scale of the baryon acoustic oscillations as imprinted in the cosmic microwave 
background and in the large scale distribution of galaxies. We have used several different 
methods to constrain the expansion history of the universe, all of which give the same 
qualitative result. 

From SN la data alone, it is evident that the universal expansion is accelerating at 
low redshifts. In particular our new sliding expansion redshift method is able to detect 
acceleration today at > 12a, even allowing for curvature. Although there are hints from 
SN la data that the universal expansion may have decelerated at high redshifts - as 
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expected in the concordance cosmological model - we can only draw that conclusion with 
high confidence if we also include CMB and BAO data, together with the assumption 
of a flat or open universe. If the universe has a small positive curvature (fifc < 0), it 
is possible, although it requires a certain level of fine tuning, to accommodate the data 
with acceleration also at high redshifts. 
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